Astronomy & Astrophysics manuscript no. 6965 


© ESO 2008 


February 5, 2008 





Visual/infrared interferometry of Orion Trapezium stars: 
Preliminary dynamical orbit and aperture synthesis imaging of the 

^^OrionisC system 

S. Kraus\ Y. Y. Balega^, J.-P. Berger^ K.-H. Hofmanni, R. Millan-Gabet"^, J. D. Monnier^ K. Ohnaka\ E. Pedretti^ 
Th. Preibisch^ , D. Schertl^ , F. P Schloerb^, W. A. Traub^ and G. Weigelti 

' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

- Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnij Arkhyz, Zelencliuk region, Karachai-Cherkesia, 357 147, 
Russia 

^ Laboratoire d' Astrophysique de Grenoble, UMR 5571 Universite Joseph Fourier/CNRS, BP 53, F-38041 Grenoble Cedex 9, France 

^ Michelson Science Center, California Institute of Technology, Pasadena, CA 91 125, USA 

' Astronomy Department, University of Michigan, 500 Church Street, Ann Arbor, MI 48104, USA 

* Department of Astronomy, University of Massachusetts, LGRT-B 619E, 710 North Pleasant Street, Amherst, MA 01003, USA 

' Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02183, USA 

Received December 19, 2006; accepted February 7, 2007 

ABSTRACT 

Context. Located in the Orion Trapezium cluster, (9'OriC is one of the youngest and nearest high-mass stars (05-07) and known to 
be a close binary. 

Aims. By tracing its orbital motion, we aim to determine the orbit and dynamical mass of the system, yielding a characterization of 
the individual components and, ultimately, also new constraints for stellar evolution models in the high-mass regime. 
Methods. Using new multi-epoch visual and near-infrared bispectrum speckle interferometric observations obtained at the 
BTA 6 m telescope, and IOTA near-infrared long-baseline interferometry, we traced the orbital motion of the 0'OriC components 
over the interval 1997.8 to 2005.9, covering a significant arc of the orbit. Besides fitting the relative position and the flux ratio, we 
applied aperture synthesis techniques to our IOTA data to reconstruct a model-independent image of the 0'OriC binary system. 
Results. The orbital solutions suggest a highly eccentricity (e x 0.91) and short-period (P x 10.9 yrs) orbit. As the current astrometric 
data only allows rather weak constraints on the total dynamical mass, we present the two best-fit orbits. Of these two, the one implying 
a system mass of 48 Mq and a distance of 434 pc to the Trapezium cluster can be favored. When also taking the measured flux ratio 
and the derived location in the HR-diagram into account, we find good agreement for all observables, assuming a spectral type of 
05.5 for e'OriCl (M = 34.0 Mq, = 39 900 K) and 09.5 for C2 {M = 15.5 Mq, T^g = 31 900 K). Using IOTA, we also obtained 
first interferometric observations on 0'OriD, finding some evidence for a resolved structure, maybe by a faint, close companion. 
Conclusions. We find indications that the companion C2 is massive itself, which makes it likely that its contribution to the intense UV 
radiation field of the Trapezium cluster is non-negligible. Furthermore, the high eccentricity of the preliminary orbit solution predicts 
a very small physical separation during periastron passage (~ 1.5 AU, next passage around 2007.5), suggesting strong wind-wind 
interaction between the two O stars. 

Key words, stars: formation - stars: pre-main sequence - stars: fundamental parameters - stars: individual: 6'OriC, (S/'OriD - 
binaries: close - techniques: interferometric - methods: data analysis 



1. Introduction 

Stellar mass is the most fundamental parameter, determining, 
together with the chemical composition and the angular mo- 
mentum, the entire evolution of a given star. Stellar evolution- 
ary models connect these fundamental parameters with more 
easily accessible, but also highly uncertain observables such as 
the luminosity and the stellar temperature. Particularly towards 
the pre-main-sequence (PMS) phase and towards the extreme 
stellar masses (i.e. the low- and high-mass domain), the exist- 
ing stellar evolutionary models are still highly uncertain and 
require further empirical verification through direct and unbi- 
ased mass estimates, such as those provided by the dynamical 
masses accessible in binary systems. Recently, several studies 
were able to provide dynamical masses for low-mass PMS stars 
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(e.g. iTamazian et al] 120021: ISchaefer et alj [20031: iBoden et al.l 
l2005h . while direct mass measurements for young O-type stars 
are still lacking. 

Furthermore, in contrast to the birth of low-mass stars, the 
formation mechanism of high-mass stars is still poorly un- 
derstood. In particular, the remarkably high binary frequency 
measured for young high-mass stars might indicate that the 
way high-mass stars are born differs significantly from the 
mass accretion scenario via circumstellar disks, which is well- 
established for low- and intermediate-mass stars. For instance, 
studies conducted at the nearest high-mass star-forming region, 
the Orion Ne bular Cluster (ONC, at a distance of 440 + 34 pc, 
I.Teffriesll2007h . revealed Lj. companions per primary (for high- 
mass stars M > IOM0, [Preibisch et al. 1999), which is signifi- 
cantly higher than the mean number of companions for interme- 
diate and low-mass stars. 
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In the very center of the ONC, four OB stars form the Orion 
Trapezium; three of which (0'Ori A, B, C) ar e known to be 
muhiple (Weigelt et al. 1999t iSchertl et al.ll2003b . e'OriD (alias 
HD 37023, HR 1896, Parenago 1889) has no confirmed compan- 
ion, although a pre liminary analysis of the radial velocity by 
IVitrichenkol (l2002al) suggests that it might be a spectroscopic 
binary with a period of ~ 20 or 40 days. 

A particularly intruiging young (< 1 Myr, iHillenbrandl 

Il997h high-mass star in the Trapezium cluster is 0'OriC (aUas 
410riC, HD 37022, HR1895, Parenago 1891). e'OriC is the 
brightest source within the ONC and also the main source of 
the UV radiation ionizing the proplyds and the M42 Hii re- 
gion. A close (33 mas) companion with a near-infrared flux ra- 
tio of ~ 0.3 between the primary (0^OriCl) and the secondary 
(0'Ori C2) was discov ered in 1997 using b ispectrum speckle in- 
terferometry (Weigelt et al.ll999l) . lDonati e t al. (2002) estimated 
the mass of 0'Ori C to be 44 + 5Mq, making it the most massive 
star in the cluster. The same authors give an effective tempera- 
ture of 45 000 ± 1 000 K and a steflar radius of 8.2 + 1.1 iJ©. 
ISimon-Dfaz et alj ( |2006|) estimated the mass independently us- 
ing evolutionary tracks and by performing a quantitative analy- 
sis of 0'OriC spectra and obtained Mevoi = 33Mq and Mjpec - 
45 + 16M0, respectively. Long series of optical and UV spectro- 
scopic observations revealed that the intensity and also some line 
profiles vary in a strictly periodic wa y. With 15.422 + 0. 002 days, 
the shortest period was reported bv IStahl et alj ( Il993h . Several 
authors interpret th is periodicity , which in the meantime was also 
detected in X-ray dGagne et al.lll997l) . within an oblique mag- 
netic rotat or model, identifyi ng 15.422 d with the rotation period 
of the star IStahl et al.1 (119961) detected a steady increase in radial 
velocity, confirmed bv iDonati et alj in 2002, which suggests a 
spectroscopic binary with an orbital period of at least 8 years. 
IVitrichenkol ( l2002bl) searched for long-term periodicity in the 
radial velocity and reported two additional periods of 66 days 
and 120 years, which he interpreted as the presence of, in total, 
three components in the system. 

Given the unknown orbit of the speckle companion, it still 
must be determined which one of these periods corresponds to 
the orbital motion of C2. Since the discovery of C2 in 1997, three 
measurements performed with bispectrum speckle interferome- 
try showed that the companion indeed undertakes orbital motion 
(ISchertl et al. l l2003l) . reaching the largest separation of the two 
components in autumn 1999 with 43+2 mas. In order to fol- 
low the orbital motion, we monitored the system using infrared 
and visual bispectrum speckle interferometry and in 2005, for 
the first time, also using infrared long-baseline interferometry. 

An interesting aspect of th e dynamical history of the ONC 
was presented by iTanI (|2004|) . He proposed that the Becklin- 
Neugebauer (BN) object, which is located 45" northwest of the 
Trapezium stars, might be a runaway B star ejected from the 
0'OriC multiple system approximately 4 000 yrs ago. This sce- 
nario is based on proper motion measurements, which show that 
BN and 0'OriC recoil roughly in opposite directions, and by 
the detection of X-ray emission potentially tracing a wind bow 
shoclfl. Three-body interaction is a crucial part of this interpre- 
tation, and C2 is currently the only candidate which could have 
been involved. Therefore, a high-precision orbit measurement 
of C2 might offer the unique possibility to recover the dynami- 
cal details of this recen t stellar ejection. However, another study 
dRodrfguez et alj|2005l) also aimed to identify the multiple sys- 



Table 2. IOTA Calibrator Stars Information 



Star 


V 


H 


Spectral 


Adopted UD diameter 








Type 


[mas] 


HD 14129 


5.5 


3.3 


G8.51II 


1.01 ±0.01" 


HD 20791 


5.7 


3.5 


G8.5III 


o.89±o.or 


HD 28322 


6.2 


3.9 


G9II1 


0.82 ±0.01" 


HD 34137 


7.2 


4.4 


K2III 


0.80 ±0.01" 


HD 36134 


5.8 


3.2 


KlIII 


1.16 ±0.02" 


HD 50281 


6.6 


4.3 


K3V 


0.77 ±0.10* 


HD 63838 


6.4 


3.6 


K2III 


0.95 ±0.01" 


HD 74794 


5.7 


3.5 


KOIII 


1.07 ±0.01" 



Not es - ° UD diameter ta ken from the CHARM2 
catalog jRichichi et aU2005l) . 

* UD diameter taken from getCal tool 
( |http : //mscweb . ipac . caltech . edu/gcWeb/ gcWeb ■ j sp P . 



tem from which BN was ejecte d and ident ified Source I as the 
likely progenitor system. Later. lOdmez et al . (2005) added fur- 
ther evidence to this interpretation by identifying Source n as a 
potential third member of the decayed system. 

2. Observations and Data Reduction 

2. 1 . Bispectrum speckle interferometry 

Speckle interferometric methods are powerful techniques for 
overcoming the atmospheric perturbations and for reaching the 
diffraction-limited resolution of ground-based telescopes, both 
at near-infrared and visual wavelengths. Since the discovery of 
0'OriC2 in 1997, we have monitored the system with the Big 
Telescope Alt-azimuthal (BTA) 6.0 m telescope of the Special 
Astrophysical Observatory located on Mt. Pastukhov in Russia. 
For the speckle observations at visual wavelengths, a 1 280x 1 024 
pixel CCD with a multialkali S25 intensifier photocathode was 
used. The near-infrared speckle observations were carried out 
using one 512x512 pixel quadrant of the Rockwell HAWAII ar- 
ray in our speckle camera, with pixel sizes of 13.4 mas (/-band), 
20.2 mas (//-band), and 27 mas (/T-band) on the sky. 

For the speckle observations at infrared wavelengths, we 
recorded interferograms of 0'Ori C and of the nearby unresolved 
star 0'OriD in order to compensate for the atmospheric speckle 
transfer function. The number of interferograms and the detector 
integration times (DITs) are listed in Table [1] 

The modulus of the Fourier transform of the ob- 
ject (visibility) was obtai ned w ith the speckle interfer- 
ometry method (Labeyry Il970l) . For image reconstruc- 
tion we used the bispectrum speckle interferometry method 
(IWeigeltill977llWeigelt & Wirnitzed[T98llLohmann et al.lll983L 



(IWeigeltlll^/ 71 1 Weigelt & > 
iHofmann & Weigeltill986l) 



' However, the more recent detection of X-ray variability in intensity 
and spectrum makes it unlikely that this X- ray emission really origi- 
nates in a wind bow shock, as pointed out bv lGrosso et al.l ( 120051) . 



2.2. IOTA long-baseline interferometry 

The Infrared Optical Telescope Array (IOTA) is a three- 
telescope, long-baseline interferometer located at the Fred 
Lawrence Whipple Observatory on Mount Hopki ns, Arizona, 
opera ting at visual and near-infrared wavelengths dTraub et al.l 
|2003|) . Its three 45 cm primary Cassegrain telescopes can be 
mounted on stations along an L-shaped track, reaching 15 m 
towards a southeastern and 35 m towards a northeastern di- 
rection. After passing a tip-tilt system, which compensates 
the atmospherically induced motion of the image, and path- 
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Table 1. Observation Log. 



Target 


Instrument 


Date 


Filter" 


Detector'' 


No. Interferograms 


Calibrators 


d 






and Configuration 


[UT] 




and Mode*^ 


Target/Calibrator 








; 

e'OriC 


BTA 6m/Speckle 


1997.784 


H 


P/DIT= 150 ms 


519/641 


e'OriD 






e'OriC 


BTA 6m/Speckle 


1998.838 


K' 


H/DIT=120 ms 


438/265 


e'OriD 






6*' Ori C 


BTA 6m/Speckle 


1999.737 


J 


H/DIT= 100 ms 


516/244 


e'OriD 






0'OriC 


BTA 6m/Speckle 


1999.8189 


G' 


S/DIT=5 ms 


500/- 


— 






e'OriC 


BTA 6m/Speckle 


2000.8734 


y 


S/DIT=5 ms 


1 000/- 


— 






e'OriC 


BTA 6m/Speckle 


2001.184 


J 


H/DIT=80ms 


684/1 523 


e'OriD 






e'OriC 


BTA 6m/Speckle 


2003.8 


J 


H/DIT=160 ms 


312/424 


e'OriD 






e'OriC 


BTA 6m/Speckle 


2003.9254 


v 


S/DIT=2.5 ms 


1 500/- 


— 






e'OriC 


BTA 6m/Speckle 


2003.928 


v 


S/DIT=2.5 ms 


2 000/- 


— 






e'oric 


BTA 6m/Speckle 


2004.8216 


v 


S/DIT=5 ms 


2000/- 


— 






6'OriC 


BTA 6m/Speckle 


2006.8 


V , R' 




~ 








e'OriC 


IOTA A35-B15-C0 


2005 Dec 04 


H 


1L7R, 2L7R 


1 1 400/8 050 


HD 14129, 


HD 36134, 


HD 34137, 














HD 50281, 


HD 63838 




0'OriC 


IOTA A35-B15-C10 


2005 Dec 02 


H 


2L7R, 4L7R 


A A r\f\ 1 A r\cf\ 

4 400/4 950 


HD 34137, 


HD 50281, 


HD 63838 


e'OriC 


IOTA A35-B15-C10 


2005 Dec 03 


H 


2L7R 


4 600/2 450 


HD 28322 






e'oric 


IOTA A35-B15-C15 


2005 Dec 01 


H 


2L7R, 4L7R 


7 250/5 000 


HD 20791, 


HD 34137, 


HD 36134 


e'OriC 


IOTA A25-B15-C0 


2005 Dec 06 


H 


2L7R, 4L7R 


5 250/4 875 


HD 28322, 


HD 34137, 


HD 36134, 














HD 74794 






e'OriD 


IOTA A35-B15-C0 


2005 Dec 04 


H 


2L7R 


800/2 800 


HD 14129, HD 36134, HD 34137, 














HD 50281, 


HD 63838 




e'OriD 


IOTA A25-B15-C0 


2005 Dec 06 


H 


2L7R, 4L7R 


1 800/4 875 


HD 28322, 


HD 34137, 


HD 36134, 



HD 74794 



Notes - Filter central wavelength and bandwidth, in nm {AJhA) - V: 545/30; G': 610/20; R': 800/60; J: 1 239/138; H: 1 613/304; 
if': 2115/214. 

*' P: PICNIC detector, H: HAWAII array, S: Multialkali S25 intensifier photocathode 

For the IOTA measurements, we used different detector read mo des to adapt to the ch anging atmospheric conditions. The two numbers in 
the given 4-digit code denote the value of the loop and read parameter i jPedretti et al.lbooil ') of the PICNIC camera. Since data taken in different 
readout modes is calibrated independently, the scattering between the data sets also resembles the typical calibration errors (see Figure|4ll. 

The dash symbol in the calibrator column indicates speckle measurements for which no calibrator was observed. 



compensating delay lines, the three beams are fed into fibers 
and coupl ed pairwise onto th e IONIC3 integrated optics beam 
combiner dBerger et al.ll2003h . The interferograms are recorded 
by temporal modulation around zero optical path delay (OPD). 
Durin g data acquisition, a fringe tracker software jPedretti et alj 
l2005h continuously compensates potential OPD drifts. This al- 
lows us to measure the three interferograms nearly simultane- 
ously within the atmospheric coherence time, preserving the 
valuable closure phase (CP) information. 

For our IOTA observations, we used four different array con- 
figurations (see Table [U, obtaining the Mv-coverage shown in 
Figure[T] 0'Ori D was observed on two different array configura- 
tions, as shown in Figure|2] During each night, we systematically 
alternated between the target star and calibrators in order to de- 
termine the transfer function of the instrument. For more details 
about the calibrator stars and the number of recorded Michelson 
interferograms, refer to Table|2] 

In order to extract visibilities and CPs from the recorded 
IOTA interferograms, we used the IDRSQ data reduction soft- 
ware. Basic principles of the algorithms i mplemented in thi s 
software package were already presented in lKraus et al. I (l2005l) . 
although several details have been refined to obtain optimal re- 
sults for fainter sources as well, such as those observed in this 
study. 

To estimate the fringe amplitude (visibility squared, V^), we 
compute the continuous wavelet transform (CWT), which de- 

^ The IDRS data reduction software can be obtained from 
[ftp : //ftp ■ mpifr-bonn .mpg . de/outgoing/skraus/idrs/ i 



composes the signal into the OPD-scale domain, providing scale 
(frequency) resolution while preserving the information about 
the fringe OPD. The turbulent Earth atmosphere introduces fast- 
changing OPD variations between the combined telescopes (also 
known as atmospheric piston), degenerating the recorded in- 
terferograms. By measuring the extension of the fringe packet 
in the CWT along both the scale and the OPD axis, we iden- 
tify the scans which are most affected by this effect and reject 
them from further processing. For the remaining s cans, we apply 
a meth od similar to the procedure presented by iKervella et al.l 
(l2004h . First, the fringe peak is localized in the CWT. In order to 
minimize noise contributions, a small window around the fringe 
peak position is cut out. Then we integrate along the OPD axis, 
yielding a power spectrum. After recentering the fringe peak po- 
sition for each scan (to compensate frequency changes induced 
by atmospheric piston), we average the power spectra for all 
scans within a dataset. In the resulting averaged power spectrum, 
we fit and remove the background contributions and integrate 
over the fringe power to obtain an estimate for V^. 

Another refinement in our software concerns the CP estima- 
tion. We found that the best signal-to-noise ratio (SNR) can be 
achieved by averaging the bispectra from all scans. The bispec- 
trum is given by the triple pro duct of the Fourier transform of the 
scans at the three baselines (iHofmann & Weigelill993b . Then, 
we use the triple amplitude to select the bispectrum elements 
with the highest SNR and average the triple phases of these ele- 
ments in the complex plane to obtain the average closure phase. 
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1999.7, BTA, 7-band 2004.8, BTA, V'-band 2005.9, IOTA, //-band 




Fig. 3. Left, Middle: Bispectrum speckle J and V'-band images reconstructed from BTA-data. Right: Aperture synthesis image of 
0'OriC reconstructed from our IOTA //-band data. Besides the false-color representation in the upper row, we show the images 
below as contours with the best-fitted positions marked with star symbols (see Table |3]l. In the image in the lower right corner, the 
restoring beam for the IOTA aperture synthesis image is shown as an inset. Over the six year interval covered by the images, orbital 
motion is clearly conceivable (1999/2004/2005: p=43/24/ 1 4 mas; 0=21471917163°). 



3. Aperture synthesis imaging 

Interpreting optical long-baseline interferometric data often re- 
quires a priori knowledge about the expected source brightness 
distribution. This knowledge is used to choose an astrophysically 
motivated model whose parameters are fitted to the measured in- 
terferometric observables (as applied in Sect.|4|. 

However, the measurement of CPs allows a much more intu- 
itive approach; namely, the direct reconstruction of an aperture 
synthesis image. Due to the rather small number of telescopes 
combined in the current generation of optical interferometric 
arrays, direct image reconstruction is limited to objects with a 
rather simple source geometry; in particular, m ultiple systems 
(for images reconstruct ed from IOTA data, see iMonnier et alj 
l2004tlKraus etaDl2005h . 

Using o ur software based on the Building Block Mapping al- 
gorithm (Ho fmann & Weigelill993h . we reconstructed an aper- 
ture synthesis image of the ^'OriC system from the data col- 
lected during our IOTA run. Starting from an initial single delta- 
function, this algorithm builds up the image iteratively by adding 
components in order to minimize the least-square distance be- 
tween the measured bispectrum and the bispectrum of the re- 
constructed image. 

The resulting image is shown in Figure [3] and provides a 
model-independent representation of our data. By combining 
the data collected during six nights, we make the reasonable as- 
sumption that the orbital motion over this interval is negligible. 



The clean beam, which we used for convolution to obtain the 
final image, is rather elliptical (see inset in Figure O, represent- 
ing the asymmetries in the Mv-coverage. 

4. IVIodel fitting 

4. 1 . Binary model fitting for 9^ Or! C 

Although the aperture synthesis image presented in the last 
section might also be used to extract parameters like binary 
separation, orientation, and intensity ratio of the components 
(Ici/Ici = 0.26, p - 14.1 mas, © = 162°), more precise values, 
including error estimates, can be obtained by fitting the mea- 
sured visibilities and CPs to an analytical binary model. 

The applied model is based on equations 7-12 presented 
in iKraus et aL I (I2005h and uses the least-square Levenberg- 
Marquardt method to determine the best-fit binary separation 
vector and intensity ratio. In order to avoid potential local min- 
ima, we vary the initial values for the least-square fit on a grid, 
searching for the global minimum. 

Since the apparent stellar diameter of fl'OriC is expected to 
be only ~ 0.2 mas at the distance of Orion, for our fits we as- 
sume that both stellar components appear practically unresolved 
to the IOTA baselines. Furthermore, we assume that the relative 
position of the components did not change significantly over the 
6 nights of observation. 

Figure |4] shows the measured IOTA visibilities and CPs 
and the observables corresponding to our best-fit binary model 
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(O-C) Orbit #1 


(O-C) Orbit #2 


Telescope 


Date 


Filter 


Flux ratio 


0° 


P 
r 


Ref. 




P 




P 








1 1 


[°] 


[mas] 




L 1 


[mas] 


L 1 


[mas] 


BTA 6m/Speckle 


1997.784 


H 


0.26 ± 0.02 


226.0 ± 3 


33 ±2 


b 


+3.0 


+0.0 


+3.0 


+0.5 


BTA 6m/Speckle 


1998.838 


K' 


0.32 ± 0.03 


222.0 ± 5 


37 ±4 


b 


+3.8 


-2.6 


+3.8 


-2.5 


BTA 6m/Speckle 


1999.737 


J 


0.31 ±0.02 


214.0 ±2 


43 ± 1 


c 


-0.9 


+ 1.5 


-0.9 


+ 1.5 


BTA 6m/Speckle 


1999.8189 


G' 


0.35 + 0.04 


213.5 + 2 


42 ± 1 


- 


-1.1 


+0.5 


-1.1 


+0.5 


D Irt Ulll/opcCIViC 


?f)nn 87^A 

ZUUU.O / JH- 


V 


U.JJ 3: U.UJ 


9 1 n n + 9 

L\J.\J III A 






-0.8 


-0.8 


-0.9 


-0.8 


BTA 6m/Speckle 


2001.184 


J 


0.29 + 0.02 


208.0 ± 2 


38 ± 1 


c 


-1.6 


-2.1 


-1.7 


-2.1 


BTA 6m/Speckle 


2003.8 


J 


0.30 ± 0.02 


199.3 ± 2 


29 + 2 




+3.9 


+0.5 


+2.8 


+0.5 


BTA 6m/Speckle 


2003.9254 


v 




199.0 ± 2 


29 ±2 




+3.6 


+ 1.3 


+3.4 


+ 1.3 


BTA 6m/Speckle 


2003.928 


v 




199.1 ±2 


29 ±2 




+3.8 


+ 1.3 


+3.6 


+ 1.3 


BTA 6m/Speckle 


2004.8216 


v 


0.34 + 0.04 


190.5+4 


24 ±4 




+4.2 


+2.4 


+4.0 


+2.4 


IOTA 


2005.92055 


H 


0.28 + 0.03 


162.74 ±2 


13.55 ±0.5 




-0.5 


+0.0 


-1.0 


+0.0 


BTA 6m/Speckle 


2006.8 


V',R 






< 15 













Notes - "* Following the conve nt ion, we measure th e position angle (PA) from north to east. 
References - IWeigeh et alJ 19991. " ISchertl et al.ll2003l 
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Fig. 1. Mv-plane coverage obtained on O^Ori C with the four 
IOTA array configurations used. The relatively strong asymme- 
try in the Mv-plane coverage mainly results from the position of 
0'OriC relatively close to the celestrial equator. 



ixlJAof = 1.35,;tc/'/dof = 1-48). The separation p, PA 0, and 
intensity ratio of this binary model are given in Table [3] together 
with the positions derived from the speckle observations. To il- 
lustrate more clearly that the measured IOTA visibilities resem- 
ble a binary signature, in Figure |5] we show a projection of the 
sampled two-dimensional Fourier plane along the binary vec- 
tor, revealing the cosine modulation corresponding to the Fourier 
transform of a binary brightness distribution. 

For the speckle data (providing a complete Fourier sampling 
up to the spatial frequency corresponding to the diameter of the 
telescope primary mirror), we determine the binary parameters 
by fitting a two-dimensional cosine function directly to the 2-D 
speckle interferogram power spectrum. 




Fig. 2. Mv-plane coverage obtained on fl^OriD with two IOTA 
array configurations. Symbols and colors as in Figure[T] 



4.2. Resolved structure around 0' Ori D: 
Potential detection of a companion 

Besides the main target of our observational programme, 
O^OiiC, during the two nights with the best seeing conditions, 
we also recorded four datasets on fl'OriD. Despite lower flux 
(e'Ori D: H=5.9, fl'Ori C: i/=4.6), the quality of the derived vis- 
ibilities and CPs seems reliable, although slightly larger errors 
must be assumed. 0' Ori D appears resolved in our measurements 
with a significant non-zero CPs signal (~ 10°) on the A35-B15- 
CO baseline. This CP indicates deviations from point-symmetry, 
as expected for a binary star We applied the binary model fit 
described in Sect. 14.11 and found the binary system with an in- 
tensity ratio of 0.14, p = 18.4 mas, and = 41° (Figure|6]l to be 
the best-fit model Or^/dof = 1.36). 

However, considering the Mv-coverage of the existing 
dataset, this solution is likely not unique, and it can not be ruled 
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Fig. 5. Projection of the sampled Fourier plane along the fitted 
binary PA of 162.74° {x - u cos 0-v sin 0), clearly revealing the 
binary signature. The dashed line shows the theoretical cosine 
visibility profile for a binary star with separation 13.55 mas and 
intensity ratio 0.28. 

Table 4. Preliminary orbital solutions, dynamical parallaxes, and 
system masses for 0'Ori C. 





Orbit #1 


Orbit #2 


p 


[yrs] 


10.98 


10.85 


T 




1996.52 


1996.64 


e 




0.909 


0.925 


a 


[mas] 


41.3 


45.0 


i 


n 


105.2 


103.7 


Q. 


n 


56.5 


56.9 


w 


[°] 


65.7 


68.2 




1.61 


1.59 




[mas] 


2.304 ± 0.066 


2.585 ± 0.074 




[pc] 


434+ 12 


387 + 11 


(Mcl + Mc2)"' 


[Mq] 


47.8 ±4.2 


44.8 ± 3.9 



Notes - The errors on the dynamical parallaxes and correspond- 
ing distances were estimated by varying the measured binary flux ratio 
within the observational uncertainties, the assumed spectral types for 
the bolometric correction by one sub-class, the extinction by ±0.2 mag- 
nitudes, and by using three different MLRs (bv lBaize & Romani|[l94a : 
iHeintzll 19781 : Id ermrcan & Kahramanlll99lh . However, the given errors 
do not reflect the uncertainties on the orbital elements a and P. Due to 
the presence of the multiple orbital solutions, it is currently not possible 
to quantify these errors reliably. 



out that other geometries, such as for inclined circumstellar disk 
geometries with pronounced emission from the rim at the dust 
sublimation radius (see e.g. iMonnier et al. I l2006h . might also 
produce the asymmetry required to fit the data. 
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Fig. 6. Visibilities and Closure Phases derived from the recorded 
IOTA interferograms for fl^OriD (points with error bars). The 
solid line shows the binary model fit, described in more detail in 
Section[ 



5. Results 

5.1. Preliminary physicai orbit of the 6^ Ori C binary system 

Our multi-epoch position measurements of the 0'OriC system 
can be used to derive a preliminary dynamical orbit. To find 
orbital solutions, we used the method described by iDocobd 
(Il985h . This method generates a class of Keplerian orbits pass- 
ing through three base points. From this class of possible solu- 
tions, those orbits are selected which best agree with the mea- 
sured positions, where we use the error bars of the individual 
measurements as weight. In order to avoid over-weighting the 



orbit points which were sampled with several measurements at 
similar epochs (two measurements in 1999.7-1999.8 and three 
measurements in 2003.8-2003.9), we treated each of these clus- 
ters as single measurements. 

In Table |4] we give the orbital elements corresponding to the 
two best orbital solutions found. As the ;if~/dof values of the 
two presented orbits are practically identical, the existing as- 
trometric data does not allow us to distinguish between these 
solutions. These orbits and the corresponding O-C vectors are 
shown in Figure [7] (see Table [3] for a list of the O-C values). 
As the ephemerids in Table |5] and also the position predictions 
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Table 5. Ephemerides for the C orbital solutions presented 
in Table H 
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Fig. 7. Top: Two best-fit apparent orbits of the ^^'OriC binary 
system (solid and dashed line). Each position measurement (see 
Table[3]) is connected to its predicted position with an 0-C line. 
Besides our orbital solutions (see Table |4|, the line of nodes 
and the predicted positions for the interval 2007.0 to 2008.5 
(increments of 0.25 yrs) are shown. North is up and east is to 
the left. Bottom: Radial velocity variations of both components 
co mputed fo r our o rbital solutions using the method presented 
by iPourbaixl (Il998l) . For the mass ratio Mci/Mci, we assume 
0.47, as obtained from our modeling of the flux ratio presented 
in Sect. 15.41 When computing the radial velocity curve, the (un- 
known) system velocity Vq represents a free-parameter corre- 
sponding to a velocity off'set, which we set to -5 kms"' for this 
plot. The radial velocity measureme nts (points with errorbars) 
were taken from lVitrichenkol ( l2002bb . 



(dots) in Figure |2] show, future high-accuracy long-baseline in- 
terferometric measurements are needed to distinguish between 
these orbital solutions. 

Potentially, additional constraints on the S'OriC binary or- 
bit could be provided by radial velocity measurements, such as 
those published by Vitrichenko ( 2002b) and in the references 
therein. However, the complexity of the 0' Ori C spectrum - in- 
cluding the line variability corresponding to the magnetically 
confined wind-shock region expected towards 0'OriC - makes 
both the measurement and the interpretation of radial velocities 
for 0'OriC very challenging. Since it is unclear whether these 
velocities really correspond to the orbital motion of the binary 
system or perhaps to variations in the stellar wind from 0' Ori C, 



Epoch 


Orbit #1 


Orbit #2 






P 


Vy 


P 




L J 


[mas] 


L 1 


[mas] 


2007.0 


100.9 


8.7 


101.8 


8.7 


2007.5 


28.6 


1.9 


-22.4 


0.8 


2008.0 


229.1 


21.7 


227.9 


23.0 


2008.5 


224.6 


30.1 


223.8 


31.0 


2009.0 


221.7 


35.1 


221.1 


35.8 


2010.0 


217.5 


40.2 


217.0 


40.5 


2011.0 


213.9 


41.6 


213.4 


41.6 


2012.0 


210.3 


40.5 


209.8 


40.2 


2013.0 


206.3 


37.5 


205.8 


37.0 


2014.0 


2014 


33.0 


200.8 


32.3 


2015.0 


194.6 


27.1 


193.7 


26.3 


2016.0 


1834 


20.1 


181.7 


19.2 


2017.0 


159.7 


12.8 


155.1 


12.0 



we did not include these velocity measurements as a tough con- 
straint in the final orbital fit, but show them together with the 
radial velocities corresponding to our best-fit orbit solutions in 
Figure |7] 

Both orbital solutions suggest that during periastron pas- 
sage, the physical separation between CI and C2 decreases to 
~ 1 .5 AU, corresponding to just ~ 30 stellar radii. Besides the 
strong dynamical friction at work during such a close passage, 
strong wind-wind interaction can also be expected. 

It is worth mentioning that besides the presented best-fit 
orbital solutions, a large number of solutions with longer or- 
bital periods exist, which are also fairly consistent with the 
astrometric measurements. However, since these orbits have 
slightly higher ;^'^/dof values than the solutions presented 
above and also correspond to physically unreasonable masses 
(Mci + Mc2 i 20 or > 130 Mq, assuming ^f=440 pc), we 
rejected these formal solutions. 

5.2. Dynamical masses and parallaxes 

Kepler's third law relates the major axis a and the orbital period 
P with the product of the system mass and the cube of the par- 
allax; i.e. (Mci + Ma) ■ - cr' IP~ (where a and n are given in 
mas, P in years, and M in solar masses). 

In order to separate the system mass and the parallax 
in absence of spectrosc opic orbital elements, the method by 
iBaize & Romanil (1 19461) can be applied. This method assumes 
that the component masses follow a mass-luminosity rela- 
tion (MLR), which, together with a bolometric correction and 
extinction-corrected magnitudes, allows one to solve for the sys- 
tem mass Mci + Mc2 and the dynamical paralla x TTdyn- When 
using the MLR bv lDemircan & Kahramanid 19911). the bolomet - 
ric correction for 05.5 and 09.5 stars by iMartins et al.l ( l2005h . 
and the extinction corrected magnitudes given in Table|6] we de- 
rive the dynamical masses and parallaxes given in Table|4] When 
comparing the distances corresponding to the dynamical paral- 
laxes derived for Orbit #1 (ddyn = ^^y^ - 434 pc) and Orbit 
#2 (t/dyn - 387 pc) with dist ance estimate s from the literature 
(e.g. £/ = 440 + 34 pc from iJeflriesI 120071: see also references 
herein), orbit solution #1 appears much more likely. The dynam- 
ical system mass corresponding to Orbit #1 is 47.8 Mq, which 
must be scaled by a factor (cZ/^fdyn)^ when distances other than 
'^dyn - 434 pc are assumed. 
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5.3. The orbital parameters in tiie context of reported 
periodicities 

Several studies have already reported the detection of periodic- 
ity in the amplitude, width, or velocity of spectral lines around 
0'OriC. This makes it interesting to compare whether one of 
those periods can be attributed to the presence of companion C2: 

P X 15.422 ± 0.002 d: By far, the best-estabHshed period- 
icity towards 0'OriC was detected in hydrogen recom- 
binati on lines and vari o us ph o tospheric and stellar-wind 
lines ('Stahl et al.' '19931 Il996t IWalborn & Nicholsl Il994t 
lOudmaijer et al. 1997). L ater, the same per iod was also 
found in the X-ray flux ( iGagne et alj Il997h and even i n 
modulations in the Stokes parameters (IWade et al.l l2006h . 
Although possible associations with a hypoth etical low- 
mass stellar companion were initially discussed (IStahl et al.l 
Il996l) . this period is, in the context of the magnetic rotator 
model, most often associated with the stellar rotation period. 
We can rule out that C2 is associated with this periodicity, as 
we do not see significant motion of C2 within the seven days 
covered by the IOTA measurements^ 

60 d < < 2 yrs, P x 120 yrs: IVitrichenkol (l2002bl) fitted 
radial velocity variations assuming the presence of two com- 
panions and determined possible periods of 729. 6/L days 
(with L an integer < 13) for the first and ~ 120 yrs for the 
second companion. Since our orbital solutions do not match 
any of these periods, we consider an association of 0'Ori C2 
very unlikely. 

f ^ 8 yrs: IStahll (1 19981) rep orted a steady increase in radial 
velocitv. lOonati et aL (1200 2;) confirmed this trend and esti- 
mated that this increase might correspond to the orbital mo- 
tion of a companion with a period between 8 yrs (for a highly 
eccentric orbit) and 16 yrs (for a circular orbit). With the 
found period of ~ 11 yrs, it is indeed very tempting to asso- 
ciate 0'OriC2 with this potential spectroscopic companion. 
However, as noted in Sect. 15.11 the set of available spectro- 
scopic radial velocity measurements seems rather inhomoge- 
neous and fragmentary and might contain observational bi- 
ases due to the superposed shorter-period spectroscopic line 
variations, as noted above. 

5.4. Nature of the 9^ Oh C components 

Most studies which can be found in literature attributed the 
whole stellar flux of 0'OriC to a single component and 
determined a wide ran ge of spectral types including 05.5 
(Gagne et al."2005'), 06 (Levato & Abll ll976l:IS im6n-Df i^t all 
[2006) . 07 (van Altenaetal. 1988), to 09 (iTrumpler 1931 . In 
order to resolve this uncertainty, it might be of importance to 
take the presence of 6''OriC2 into account. Besides the con- 
straints on the dynamical masses derived in Section 15.11 addi- 
tional information about the spectral types of ff'OriCl and C2 
can be derived from the flux ratio of the components in the 
recorded bands. 

In co ntrast to our earlier studies ( IWeigelt et al.l 119991 : 
ISchertl e t al. 2003), we can now also include the V-band flux ra- 
tio measurement to constrain the spectral types of the individual 
components. The V-band is of particular interest, as a relative 
increase of the flux ratio FqiIFqi from the visual to the near- 
infrared would indicate the presence of circumstellar material, 
either as near-infrared excess emission or intrinsic extinction to- 
wards C2 (assuming similar effective temperatures for both com- 
ponents). Our speckle measurements indicate that FqiIFcx stays 



0.5 



0.4 



0.3 



0.2 - 



0.1 



I:: 



Measurements ^cd^c-t 
06/09.5 - Mo, =31 M_„, Mc2=16 M_„ (M05) 
05.5/09.5 - Mo, =34 M,,,„, Mf.2=16 M,,,„ (M05) 
05/09.5 - Mo =38 m'"" Mo2=16 m'"" (M05) 
05.5/09 - Mo, =40 M,,,„, Mo2=16 M,,,„ (C04) 
07/B1-Mg;=32 M^:, Mg=12M^:ic04i 



0.4 



0.6 



1.2 1.4 1.6 

Wavelength [^m] 



2.2 



Fig. 8. Measured intensity ratio of the 0'OriC components as 
a function of wavelength (points with errorbars). For various 
spectral-type combinations, the curves show the expected in- 
tensity ratio, assuming the stars contribute purely photospheric 
emission (black-body emission with luminosities and effec- 
tive temperatures as given in t he stellar evol ution models from 
iMartins et al.l(l2005L M05) and lciarei (|2004 C04)). 



rather constant from the visual to the near-infrared. Therefore, in 
the following we assume that the major contribution of fl'Ori C2 
to the measured flux is photospheric. 

In Figure |8] we show the measured Fci/Fci as a function 
of wavelength and compare it to model curves corresponding to 
various spectral-type combinations for CI and C2. To compute 
the model flux ratios, we simulate the stellar photospheric emis- 
sion as black-body emission B{T f.ff) with eff'ective temperatures 
Teif and stellar radii R, as predicted b y stellar evolutionary mod- 
els ( IClarelll2004HMartins et al.ll2005l) : 



(1) 



Under these assumptions, the companion C2 would have to be 
rather massive {MqiIMqi - 0.45 + 0.15) to obtain reasonable 
agreement with the measured flux ratios (see Figure|8]l. 

Using a value for Ay from Uterature, the flux ratios can also 
be used to estimate the photometry of the individual components 
(Table |6]l. Then, the spectral type of CI and C2 can be deter- 
mined by comparing the location of the stars in the HR-diagram 
with stellar evolution mo dels. For this, we adopt the procedure 
from ISchertl et al.l (l2003l) and convert the derived photometry 
into locations in the HR-diagram usin g the colors and bolomet- 
ric correcti ons from | Kenvon & Hartman n ( 1995, and references 
therein) and lMartins &Plezl(l2006ir Assuming coevality for both 
stars, the spectral type of the individual components can be con- 
strained by finding the location where the curves for the vari- 
ous spectral bands and the isochrone intersect. As can be seen 
in Figure |9] the allowed locations for CI intersect the Zero- 
Age Main SequenctO (ZAMS) around Jeff = 46 000 + 4 000 K, 
logL/LQ = 5.3 + 0.2 (coiTesponding to 05) and around Teff = 
33 000 + 2 000 K, log L/L© = 4.5 + 0.1 (corresponding to 09) 
for C2. 



^ With a dynamical age of ~ 3 x 1 0^ yrs, it seems ju stified that the 
Trapezium stars are real ZAMS stars dSchulz et alj|2003l) . although the 
strong magnetic activi ty from 0'OriC w as also associated with a pre- 
main-sequence origin jPonati et alj200A . 
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Table 6. Derived d ereddened magnitudes a nd colors for the 0'OriC compo nents. For the pho tometry of the total system, 
we used data from 'Hillenbrand^etalJ (1 19981 7=4.63, //=4.48, K^4A1) and 'Hillejibrand| (Il997l y=5.12). A n extinction of 



Av - 1.74 dHille nbrand, 1997.) was assumed (using the reddening law by iMathis., 1 990i and, similar to 
high 7? V = 5.0) 



Mathis & Wallenhorst,198L a 
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Fig. 9. HR-diagram with PMS evolutionary tracks (Z = 0.02; 
thin solid lines with the corresponding masses lab eled aside). 
For masses < 7M(T| (> 9M(7)), the tracks from ISiess et all 
I2OOOI (iBernasconi & Maededll996h were used. The three dot- 
ted lines represent the isoch rones for 0.3, 0.6, and 1.0 Myr 
(IBernasconi & Maedej|1996l) . Usi ng the synthetic colors and 
bolometric correctio ns compiled by Kenvon & Hartrnannl(ll995l 
including data from |Johns"onlll966[ iBessell & Brett 1 19881 and 
others) and th e O s tar colors by iMartins et alj ( 2005 ) and 
iMartins & Pled ( l2006h . we converted the measured photometry 
for 0' Ori CI, C2, Dl, and D2 into the corresponding allowed 
locations in the HR-diagram (thin red lines). 



We conclude that the spectral type combination, which si- 
multaneously provides good agreement to the measured flux 
ratios, the HR-diagram, and the dynamical masses derived in 
Sect. 15.11 is given by the following stellar par ameters (using the 
evolutionary models from IMartins et al.1l2005h : 
CI: 05.5 (M = 34.0 Mq, T^s = 39 900 K, logL/L© = 5.41) 
C2: 09.5 (M = 15.5 M©, Teff = 31 900 K, log L/L© = 4.68) 

5.5. Nature of the potential 0^ Ori D companion 

Although the 0'OriD binary parameters presented in Sect. 14.21 
must be considered preliminary, it might be interesting to de- 
termine the spectral type of the putative components. We apply 
the procedure discussed in Sect. 15.41 to determine the photom- 
etry of the components from the measu red intensity ratio (pho- 
tometry for the unresolved system from lHillenbrand et al.lll998t 
//=5.84) and derive ///5i=5.98 and //o2=8.12, respectively. 



Searching again for the intersection between the allowed lo- 
cations in the HR-diagram with the isochrones applicable to the 
ONC (Figure O, the best agreement for Dl can be found with 
Teif = 31 500 + 4000 K, log L/L© = 4.25 + 0.1 (corresponding 
to 09. 5). Accordingly, D2 might be either a B4 or B5 type star 
which has just reached the ZAMS (T^ff = 16 000 + 4000 K, 
logL/L© - 2.6 ± 0.2) or a pre-main-sequence KO type star 
(r, ff = 5 000 ± 1 000 K log L/L© = 1.3 + 0.2). 

IVitricherikol ( l2002al) examined radial velocity variations of 
0'OriD and presented preliminary spectroscopic orbital ele- 
ments for a companion with a 20.2 d period (or twice that pe- 
riod, P=40.5 d). Assuming 20 M© as the system mass, these 
periods correspond to a major axis of 0.05 or 0.08 AU (~ 0.1 
or 0.2 mas). Since this is far below the 18 mas suggested by our 
binary model fit, we do not associate our potential companion 
with the proposed spectroscopic companion. 

The multiplicity rate in a young stellar population such as 
the Trapezium cluster is an important quantity which might al- 
low us to draw conclusions not only about the dynamical history 
of the ONC, but also about the mechanisms controlling the star 
formation process. The detection of a new companion around 
0'Ori D further increases the multiplicity rate for high-mass stars 
in the ONC. For instance, c onsidering the sample of 13 Orion O- 
and B-type stars studied bv lPreibisch et al.l (Il999h now yields 10 
visual and 5 spectroscopic detected companions (including one 
quintuple system, namely 0'OriB). This corresponds to an aver- 
age observed companion star frequency (CSF) of 1.15 compan- 
ions per primary. Despite the fact that this value only represents 
a strict lower limit due to observational incompleteness, it is al- 
ready higher than the incompleteness-corrected CSF determined 
bv lDuq uennov & Mayor ( 1991) for a distance-limi ted sample of 
solar-t ype field stars (0.5 companions per primary). iKohler et al.l 
(I2OO6I) have reported that the CSF for low- and intermediate- 
mass stars in the ONC is about a factor of 2.3 lower than the 
CSF in the lDuquennov & Mavoj sample, making the differences 
in the CSF between the low-, intermediate-, and high-mass star 
population in the ONC highly significant. Several studies (e.g. 
IPr eibisch et"ai]|1999l: iBallv & Zinnecke^l2005l: iBonnell & Bate! 
200^ have already interpreted this as evidence that diff'erent for- 
mation mechanisms (e.g. stellar coalescence vs. accretion) might 
be at work in different mass regimes. 



6. Conclusions 

We have presented new bispectrum speckle interferometric and 
infrared long-baseline interferometric observations of the Orion 
Trapezium stars 0'OriC and D. This data was used to recon- 
struct diffraction-limited NIR and visual speckle images of the 
fl'OriC binary system and, to our knowledge, the first model- 
independent, long-baseline aperture-synthesis image of a young 
star at infrared wavelengths. 

For 0'Ori D, we find some indications that the system was re- 
solved by the IOTA interferometer Although the non-zero clo- 
sure phase signal suggests asymmetries in the brightness dis- 
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tribution (maybe indicative of a close companion star), further 
observations are required to confirm this finding. 

From our multi-epoch observations on 0'OriC (covering 
the interval 1997.8 to 2005.8), we derived the relative position 
of the companions using model-fitting techniques, clearly trac- 
ing orbital motion. We presented two preliminary orbital solu- 
tions, of which one can be favoured due to theoretical argu- 
ments. This solution implies a period of 10.98 yrs, a semi-major 
axis of 41.3 mas, a total system mass of ~ 48M0, and a dis- 
tance of 434 pc. Furthermore, we find strong indications that 
0'OriC2 will undergo periastron passage in mid 2007. As the 
binary separation at periastron is expected to be ~ 1 mas, fur- 
ther long -baseline interferometric observations on fl'OriC are 
urgently needed to refine the orbital elements, the stellar masses, 
and orbital parallaxes. Through comparison with stellar evolu- 
tionary models and modeling of the measured intensity ratio, 
we find evidence that the companion 6?'OriC2 is more massive 
(Mci/Mci ~ 0.45 + 0.15) than previously thought; likely of late 
O (09/9.5) or early B-type (BO). The contribution of the com- 
panion to the total flux of S'OriC and the interaction between 
both stars might be of special importance for a deeper under- 
standing of this intriguing object. Therefore, we strongly encour- 
age observers to acquire high dispersion spectra of the system 
in order to trace the expected radial velocity variations and the 
wind-wind interaction of the system. 
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Fig. 4. Visibilities and Closure Phases derived from the recorded IOTA interferograms forfliOriC as a function of hour angle (points 
with error bars). The solid line shows the binary model fit, described in more detail in Section |4] The different symbols represent 
the different detector modes used (see Table [TJ. The data for each detector mode was calibrated separately, so the scattering of the 
data groups represents the typical caUbration errors. 



